Triaxial nuclear models and the outer crust of nonaccreting cold neutron stars 
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The properties and composition of the outer crust of nonaccreting cold neutron stars are stud- 
ied by applying the model of Baym, Pethick, and Sutherland (BPS) and taking into account for 
the first time triaxial deformations of nuclei. Two theoretical nuclear models, Hartree-Fock plus 
pairing in the BCS approximation (HF-BCS) with Skyrme SLy6 parametrization and Hartree-Fock- 
Bogolyubov (HFB) with Gogny D1S force, are used to calculate the nuclear masses. The two 
theoretical calculations are compared concerning their neutron drip line, binding energies, magic 
neutron numbers, and the sequence of nuclei in the outer crust of nonaccreting cold neutron stars, 
with special emphasis on the effect of triaxial deformations. The BPS model is extended by the 
higher-order corrections for the atomic binding, screening, exchange and zero-point energies. The 
influence of the higher-order corrections on the sequence of the outer crust is investigated. 

PACS numbers: 26.60.+C, 21.10.Dr, 21.60.-a, 97.60.Jd 



I. 



INTRODUCTION 



Neutron stars are believed to be born as the central 
compact remnant of an ordinary star which exploded 
in a core-collapse supernova. With initial temperatures 
above 10 11 K one can assume that the matter of the neu- 
tron star is at local nuclear statistical equilibrium. For 
a nonaccreting neutron star it is plausible that its com- 
position is determined by this condition even when the 
neutron star has cooled down to much lower tempera- 
ture (here we are ignoring effects from possible contami- 
nations of fall-back supernova material) . 

Within the first few centimeters of the outermost lay- 
ers of a cooled-down neutron star the density already 
becomes so high that atoms begin to touch. At p ~ 10 4 
g/cm 3 the atoms are completely ionized and separated 
into their constituents, the electrons and nuclei. The 
light electrons form a degenerate Fermi gas, soon becom- 
ing relativistic (at p ~ 10 7 g/cm 3 ), whereas the massive 
nuclei arrange in a solid (bcc)-lattice to minimize the 
Coulomb-interaction. To determine the ground state nu- 
cleus in the sequence of the outer crust of nonaccreting 
neutron stars, it is essential to know the binding ener- 
gies of all nuclei from the valley of /3-stability up to the 
neutron drip line. With the lowest energy per nucleon, 
56 Fe is present at low density. With increasing baryon 
and electron chemical potential (whereas the latter is 
fixed by charge neutrality) the ground-state nuclei be- 
come more and more neutron-rich until finally no more 
neutrons can be bound. This happens at the neutron drip 
(pnd ~ 4.3 x 10 11 g/cm 3 ), where free neutrons begin to 
appear and the inner crust begins. 

In 1971 Baym, Pethick and Sutherland (BPS) Q ana- 
lyzed this form of cold dense matter and calculated the 
resulting sequence of nuclei and the equation of state of 
the outer crust of nonaccreting cold neutron stars. They 
used the nuclear data from the droplet model of My- 
ers and Swiatecki [2| for their analysis. Over the years 
more and more masses of unstable nuclei were measured 
experimentally. In parallel the theoretical nuclear mod- 



els, needed for nuclei with unknown mass, have been 
developed further. Several publications followed, which 
updated the results of BPS by using the latest experi- 
mental and theoretical nuclear data: In 1989 Haensel, 
Zdunik and Dobaczewski [3[ used a Skyrme Hartree- 
Fock-Bogolyubov (HFB) calculation with the parameter 
set SkP in spherical approximation [J] , and a newer ver- 
sion of the droplet model from Myers Q. Haensel and 
Pichon 6] in 1994 used the experimental nuclear data 
from the 1992 atomic mass table of Audi and Wapstra 
[7| and the theoretical nuclear mass tables of the droplet 
models from Moller and Nix Q, as well as the Skyrme 
parametrization of Aboussir et al. @. The last update 
was done by Riister, Hempel and Schaffner-Bielich [T(| by 
applying the 2003 atomic mass table from Audi, Wapstra, 
and Thibault [ll| and a comprehensive set of 21 differ- 
ent theoretical nuclear models. Besides the finite-range 
droplet model (FRDM) [l2|, EH various non-relativistic 
Skyrme-parameterizations and mass tables based on rel- 
ativistic nuclear field theories were employed. Further- 
more the effects of pairing and axial deformations were 
examined. In the present work, for the first time the im- 
pact of triaxial deformations on the neutron drip line and 
the sequence of nuclei of the outer crust will be studied. 
Haensel and Pichon Q extended the original BPS model 
by including higher-order corrections of the atomic bind- 
ing, screening, and exchange energies. In the present 
calculation these corrections are also applied and their 
influence on the equation of state and the sequence of 
ground state nuclei will be analyzed. 



The paper is organized as follows. Section |TT] gives 
a brief outline of the two triaxial nuclear models. The 
BPS model and its higher-order corrections are present in 
Sec. IIIII The binding energies, the location of drip lines, 
and the sequences of the outer crust between the two 
theoretical models are compared in Sec. IIVI The effect 
of higer order corrections is also discussed. Section [V] is 
devoted to the summary and conclusion. 
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II. DESCRIPTION OF THE TRIAXIAL 
NUCLEAR MODELS 

In this section, we describe the two mean-field meth- 
ods [lH used in the present studies: a Skyrme energy 
functional with BCS pairing and a HFB theory with the 
finite-range Gogny interaction. Both calculations include 
nuclear triaxial features, where the gamma degree of free- 
dom plays an important role in the description of de- 
formed nuclei. 



A. Skyrme Hartree-Fock plus BCS 

Density functional theory states that the total energy 
of many-body system can be formulated by an energy 
functional 

£ = £kinM + £sk[p, T, J] + £c[Pp] ~ £c.m., (1) 

which only depends on the local distribution of density 
p, kinetic energy density r and spin-orbit current J. 
Here f kin is the functional of the kinetic energy, £sk the 
Skyrme energy functional, £q the Coulomb functional in- 
cluding the exchange term in Slater approximation, and 
£ c .m. the center-of-mass correction [15|, [l6(. The pairing 
correlations are treated in the BCS approximation using 
a delta pairing force [HI EH, V va i T (r,?) = V q 5(r - ?). 
The pairing strength V p for the protons and V n for the 
neutrons are fitted to the pairing gaps in isotopic and iso- 
tonic chains [l9j |. The pairing energy functional is given 
by 

fpair = \j2 V « f d " rA l ( 2 ) 
q— p,n 

with A g the pairing density. 

The variation of the energy functional £ with respect 
to the single-particle wave functions yields the mean-field 
equations 

d£ 

h<j> k = € k (j) k with h = — . (3) 

The coupled HF-BCS equations are solved on a grid 
in coordinate space with the damped gradient iteration 
method [20| and a Fourier representation of the deriva- 
tives. No symmetry restriction has been imposed in the 
calculation. 

In the present work, we have chosen the SLy6 
parametrization [2l| from recent fits. The fitting of the 
force SLy6 laid particular emphasis on the properties of 
neutron matter and neutron rich nuclei in order to im- 
prove the isospin properties away from the /3-stability 
line. Nuclear matter properties, the binding energies and 
radii of the doubly closed-shell nuclei 40 ' 48 Ca, 56 Ni, 132 Sn, 
and 208 Pb were used for the fit. The SLy6 force is a good 
candidate for describing isotopic properties of nuclei from 
the /3-stability line to the drip line. 



B. Gogny Hartree-Fock-Bogoliubov 

The other method used in the present study is HFB 
with finite-range Gogny interaction [22|, [H, [Hj]. The 
HFB equation 

(-A* -t*)(S) =£fc (S)' W 

is expressed in terms oi h = t + T — A with T being the HF 
potential and A the pairing potential. The total energy 
of the nuclear system Eq reads 

E = Tr(fp) + iTr(r P ) - ^Tr(A«;*), (5) 

with the density matrix p and pairing tensor k. 

The HFB equation was solved in a three-dimensional 
harmonic-oscillator basis [24|, 1251, |26j. The triaxial os- 
cillator parameters in the Hermite polynomials were op- 
timized for each nucleus to maximize the ground-state 
binding energy. All the contributions to the HF and 
pairing fields arising from the Gogny and Coulomb in- 
teractions as well as the two-body correction of the ki- 
netic energy are included in the self-consistent proce- 
dure. Here we employed the finite-range Gogny force 
with parametrization D1S [13, H^j, which was adjusted 
to give a better description of the nuclear surface energy. 
The finite-range Gogny force provides both the HF mean- 
field and pairing field simultaneously in the framework of 
the full HFB theory. 



III. THE EXTENDED BPS MODEL 

To calculate the outer crust of nonaccreting cold neu- 
tron stars with the two triaxial nuclear models, we start 
from the BPS model [l[ specified in [l(| ■ The total energy 
density 

etot = n N {W N + W L ) + e e , (6) 

describes completely ionized nuclei with mass Wn, and 
number density un immersed within an ideal Fermi gas 
of electrons (e e ) and the corresponding Coulomb interac- 
tion. The nuclei are arranged in a body-centered cu- 
bic lattice, represented by the lattice energy Wl- In 
the present work the energy density Eq. ([6j) is extended 
with the higher-order corrections of atomic binding (i? e /), 
screening (Ws c ), exchange [We x ), and zero-point energy 
(Wzp), 

etot = n N (W N + W L +B el + W S c + W E x + Wzp)+e e . (7) 

For the masses of the nuclei Wn, which are the essen- 
tial input information for the calculation, we always pre- 
fer taking the experimental data of the atomic mass table 
2003 from Audi, Wapstra, and Thibault [ll|, if available. 
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As nuclear binding energies are needed, the tabulated val- 
ues are corrected for the included atomic electron binding 
energy with an empiric formula given in 29] 

B el = (14.4381Z 2 39 + 1.55468 ■ 1(T 6 Z 5 - 35 ) eV , (8) 

with the charge number Z. The predicted theoretical 
binding energies do not include any atomic interaction, 
hence this correction is not necessary for them. 
The screening or Thomas-Fermi energy [30l | 



162 



4 \ 2/3 

b 175 \9ir J P 



(9) 



includes deviations of the electron distribution from uni- 
formity caused by the positively charged ions, where /i e 
is the electron chemical potential and a the fine structure 
constant. 

The third correction, the exchange energy [3C| 



where 



W Ex = -Z—am e c 2 x4>(x) , 

47T 



= k F Jr 



(10) 



(11) 



and 
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FIG. 1: Binding energies per nucleon of Ge isotopes from 
the proton drip-line to the neutron drip-line with Gogny and 
SLy6 forces. The lower panel shows the energy differences to 
FRDM values [l2l ] and the upper panel the differences to the 
experimental data [3lj . 



with j3 = x + \/l + x 2 and the mass m e and the Fermi 
momentum kp c of the electrons, takes the fermionic na- 
ture of the electron-electron interaction into account. 
The zero-point motion of the nuclei in the lattice [30l ] 
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(13) 



is much smaller than the first three corrections, and is 
only included to verify the stability of the lattice. Its 
influence on the composition and the equation of state is 
negligible. 



IV. DISCUSSION AND RESULTS 

We systematically calculate the properties of ground 
states from stable to neutron-rich nuclei with proton 
number Z = 20 to 50 using two triaxial nuclear mod- 
els, HFB with Gogny D1S parametrization and HF-BCS 
with Skyrme SLy6 parametrization. Both calculations 
predicted that most of Germanium isotopes are triaxially 
deformed, where the calculated quadrupole and triaxial 
deformations are in good agreement with the available 
experimental data [321 ]. Triaxiality is also predicted for 



some other nuclei with small differences between the two 
theoretical nuclear models. As an example, we select the 
typical triaxial deformed Ge isotopes and compare the 
calculated energies of two theoretical models with the 
available experimental data and the hnite-range droplet 
model (FRDM). Figure [1] shows the energies per nucleon 
of the Ge isotopes from the proton-rich isotope with neu- 
tron number N = 26 to neutron-rich nuclei up to the 
neutron number N = 76 which is beyond the neutron 
drip- line of all sets discussed here (see below). Owing to 
a huge energy variation along the long isotopic chain, we 
show the relative energies per nucleon, in the upper panel 
relative to the experimental data and in the lower panel 
relative to the results from the macroscopic-microscopic 
model FRDM taken from [l2j . In the range of experimen- 
tally measured isotopes, Gogny HFB calculation underes- 
timates the binding energy throughout. The deviation of 
the results for SLy6 HF-BCS from the experimental data 
is small which is no surprise because SLy6 is a more re- 
cent parametrization. The lower part of the figure shows 
the relative energies in a much broader range of isotopes 
with respect to the FRDM as experimental data is not 
available in the neutron-rich regime. The trends seen in 
the range of experimentally accessible nuclei are basically 
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FIG. 2: Two neutron separation energies for Gogny, SLy6, 
FRDM, and the available experimental data. 



continued for neutron rich nuclei. The Gogny HFB cal- 
culation tends to lower binding energies while the SLy6 
HF-BCS shows relatively small deviations in comparison 
to the FRDM. 

Figure [2] shows the two neutron separation energies 
for Gogny HFB, SLy6 HF-BCS, FRDM, and the avail- 
able experimental data. The overall agreement of the 
two-neutron separation energies is fair among different 
models and the experimental data, especially in the ex- 
perimentally measured regime up to A = 82, although 
Gogny HFB calculation systematically underestimates 
the binding energies. The location of drip line from the 
Gogny HFB calculation is different from that of SLy6 and 
FRDM, e.g. Gogny HFB predicts the location of drip- 
line at the nucleus 102 Ge while the SLy6 HF-BCS and 
FRDM predicts 106 Ge as being the last stable nucleus. 

As shown in 10], the equation of state of the outer 
crust of nonaccreting cold neutron stars is almost inde- 
pendent of the theoretical nuclear model used, and this 
remains valid for the presently examined triaxial models. 
The transitions from one equilibrium nucleus to another 
happen at different densities, while the overall form of 
the equation of state does not change. The higher-order 
corrections affect the equation of state only very slightly. 
Electronic versions of the equation of state and the chem- 
ical composition for the model SLy6 HF-BCS and some 
other selected nuclear models described in [lOj, can be 
downloaded as tables from one of the authors' webpages 
1 or the EPAPS online archive 2 . All of these tables in- 
clude the higher order corrections. 



1 http : //www. th.physik.uni-f rankf urt . de/~hempel 

2 See EPAPS Document No. E-PRVCAN- 76-004801 for 
electronic versions of the equation of state and chemi- 
cal composition in tabular form for six different selected 
nuclear models. For more information on EPAPS, see 
http: / / www.aip.org/pubservs / epaps.html 
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FIG. 3: Fractions of the energy density of the electrons e e , the 
lattice el, and the sum of the higher-order corrections ehoc on 
the total energy density e to t, calculated with the theoretical 
nuclear model SLy6 HF-BCS. 
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FIG. 4: Fractions of the pressure of the lattice Pl, and the 
sum of the higher-order corrections Phoc on the total pressure 
Ptot, calculated with the theoretical nuclear model SLy6 HF- 
BCS. 



Figure [3] shows the different contributions to the to- 
tal energy density, which is dominated by the mass of 
the nuclei. The energy density of the electrons is two 
to three orders of magnitude smaller, followed by the 
negative lattice energy density cl = WlUn, which is 
again one to two orders of magnitude below the electron 
contribution. The sum of the higher-order corrections 
e/ioc = n^{B e i + Wsc + Wex + W zp ) is positive for densi- 
ties below 10 7 g/cm 3 and negative for higher mass den- 
sities. The highest fraction of the order of 10~ 6 to 1CC D 
on the total energy density is reached at high densities 
above 10 9 g/cm 3 . The jumps in the curves arise from the 
transitions of one ground state nucleus to another. 

In Fig. 2] the contributions to the total pressure are 
depicted. The pressure is mainly generated by the elec- 
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FIG. 5: Nuclear chart (proton number Z, neutron number N) 
of the nuclei considered in this paper, taken from the atomic 
mass table [llj (crosses) and calculated results with the SLy6 
HF-BCS model (dots). The open circles mark triaxially de- 
formed nuclei. The thick line with full circles shows the se- 
quence of ground state nuclei in the outer crust of nonaccret- 
ing cold neutron stars with increasing density, the thin line 
the neutron drip line of the theoretical nuclear model. The 
dashed line with squares represents the sequence without the 
higher-order corrections. 
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FIG. 6: As Fig. \E\ now for the Gogny HFB calculation. 



trons. The lattice gives a significant contribution of 
about 40% at low densities which decreases down to 
about 4% at high densities. The pressure of the higher- 
order corrections is always negative, besides the little 
bump around 10 s g/cm 3 . The biggest influence of the 
higher-order corrections on the equation of state is seen 
at smaller densities: they lower the pressure up to ~ 3%. 
Above p ~ 10 6 g/cm 3 their contribution drops below 
1CP 3 . At densities larger than 10 9 g/cm 3 the contribu- 
tion raises again to fractions of 10~ 4 to 10 -3 . 

Figures [5] and [5] show all the nuclei examined in this 
work, using the atomic mass table and the two theoreti- 



cal models studied here. Nuclei in the theoretical models 
which fulfill the criteria 0.03 < /3 and 5° < 7 < 55° for 
the deformation parameters (3 and 7 are marked addition- 
ally with an open circle as being " triaxially deformed" . In 
both models almost all the Germanium-isotopes (Z = 32) 
are triaxially deformed. Triaxial nuclei are also found 
in the region of 56 < N < 74 and 38 < Z < 46. The 
Gogny HFB model shows only a few other triaxial nuclei, 
while for the SLy6 HF-BCS model there are some more at 
Z = 30 and in the region of Z > 44, and 96 < N < 108. 

The drip lines of the two theoretical calculations as de- 
termined via the two-neutron separation energy are also 
depicted. For both models one recognizes a rather linear 
increase at small N, followed by a pronounced vertical 
step at the magic neutron number N = 82. In compar- 
ison to the non-relativistic "state-of-the-art" theoretical 
models in [l(J, e.g. BSk8, SLy4, and FRDM, the drip 
line of the SLy6 HF-BCS model shows a similar behavior 
up to N — 82. For larger neutron numbers its drip line is 
located at AZ = — 4 larger values of proton number Z. 
The drip line of the triaxial HF-BCS calculation is also 
similar to the axial Skyrme calculations (33[, and is only 
shifted by one or two nuclei for some isotopes. If one 
compares the drip line of SLy6 HF-BCS to the Gogny 
HFB, one sees that the Gogny HFB drip line is shifted 
about AA" = — 6 to smaller neutron numbers. 

The sequence of nuclei which form the ground state 
in the outer crust of nonaccreting neutron stars with in- 
creasing density is much more sensitive to the underlying 
theoretical nuclear model than the equation of state. Fig- 
ures [5] and [B] show the results for the sequences of the two 
triaxial calculations. The last nucleus of the sequence 
with experimentally determined mass is 80 Zn which both 
models have in common. For the HF-BCS model at den- 
sities above 4.9 x 10 10 g/cm 3 the theoretically calculated 
nucleus 78 Ni with the neutron magic number 50 appears 
in the sequence. Afterwards the sequence continues along 
the pronounced neutron magic number 82 from Ru to 
fgSr, where the neutron drip line is reached. Compared 
to the previous calculations in [l(| , the sequence of the 
triaxial HF-BCS model is the same as that of FRDM 
[T2I [Hj], except that the sequence of HF-BCS is only 
shifted by AZ = 2 to nuclei with larger ZatiV = 82. 

The Gogny HFB calculation shows some distinctive 
features in comparison to the triaxial HF-BCS model and 
the theoretical nuclear models studied in After 80 Zn 
the sequence jumps to 13 gCd, as seen in Fig. which 
still originates from the atomic mass table [11| ■ None of 
the theoretical models examined in [To| predicted such a 
high proton number nucleus in the sequence. Afterwards 
a gap in the N — 82 isotone chain appears, until the first 
experimentally unknown nucleus JfRu shows up in the 
sequence. The following N = 82 sequence is the same as 
in the triaxial HF-BCS calculation. 

This behavior may be attributed to the fact that in 
the analyzed region the nuclei of the Gogny HFB cal- 
culation are more weakly bound compared to the nuclei 
with measured mass and to the other theoretical mod- 
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els. Also the location of the drip line enhances the view 
that the Gogny HFB calculation predicts relatively low 
binding energies for neutron-rich nuclei. The parameteri- 
zation D1S was constructed with emphasis on the surface 
energy and not particularly for describing binding ener- 
gies of neutron rich nuclei, which might be one of the 
reasons for the observed deviations from the other nu- 
clear sequences. But even more important could be that 
the calculation was performed in the harmonic oscillator 
basis. 

In order to show the effect of higher order correc- 
tions, Figures and [5] depict the sequences where the 
higher-order corrections were neglected. As expected, the 
higher-order corrections only lead to small changes in the 
sequence, which occur at high densities (p > 1.4 x 10 11 
g/cm 3 ). For the HF-BCS model the nucleus Ru drops 
out in the sequence, while for Gogny HFB 78 Ni enters in 
after 13 gCd. The attraction caused by the higher-order 
corrections slightly favors nuclei with larger proton and 
mass number. Without the higher-order corrections, as 
done in [Io( ]. the sequence of the HF-BCS model is more 
similar to the one of the FRDM [10]. The Gogny HFB 
calculation shows the uncommon features, namely the 
appearance of large proton numbers and the jump in the 
sequence. 

None of the nuclei appearing in the outer crust is 
triaxially (and not even axially) deformed. Deforma- 
tions become important in the regions far from closed 
shells, whereas the sequence runs mainly along the neu- 
tron magic numbers. The effect of the additional degree 
of freedom of triaxial deformation is small compared to 
the enhanced binding energy of magic nuclei, thus the se- 
quence shows no significant differences to most of the ax- 
ially deformed or even spherically calculated nuclei [Iol |. 



V. SUMMARY 

In the present work we have investigated for the first 
time the impact of triaxial deformations on the neutron 
drip line and the sequence of nuclei in the outer crust of 
nonaccreting cold neutron stars. Two theoretical nuclear 
models HF-BCS and Gogny HFB, including nuclear tri- 
axial feature, have been used to systematically study the 
ground state binding energies and triaxial deformations 
for neutron-rich nuclei from proton number Z—20 to 50. 
These neutron-rich nuclei will become experimentally ac- 
cessible in the near future with the ISAC-U and upcom- 
ing facilities as FAIR at GSI Darmstadt. Taking the 
triaxially deformed Ge isotopes as an example, we found 
that the Gogny HFB calculation systematically underes- 
timates the binding energies in comparison to the FRDM 
results and to the experimental data, while the HF-BCS 
calculation shows smaller deviations. The two-neutron 
separation energies among Gogny HFB, SLy6 HF-BCS, 
FRDM and the experimental data are in good agreement 
around the valley of stability, while the Gogny HFB cal- 
culation predicts the drip-line nucleus to be 102 Ge com- 



pared to 106 Ge for SLy6 HF-BCS and FRDM. 

The location of the drip line of the triaxial HF-BCS 
calculation is similar to that of the axial Skyrme calcu- 
lation [33j . and is shifted by one or two nuclei for some 
isotopes. In comparison to the non-relativistic theoreti- 
cal models in [T(J, e.g. BSk8, SLy4, and FRDM, the drip 
line of the HF-BCS calculation shows a similar behavior 
up to N — 82, while for larger neutron numbers it is lo- 
cated at somewhat larger values of Z. The drip line of the 
Gogny HFB calculation lies at smaller neutron numbers 
compared to the HF-BCS case, indicating overall slightly 
lower binding energies. In both models the magic neu- 
tron number N = 82 is observed as a pronounced vertical 
step in the drip lines. 

We presented the sequence of ground state nuclei for 
the outer crust of neutron stars with an extended BPS 
model. The sequence of the HF-BCS calculation is sim- 
ilar to those obtained in [l(| for spherical or axially de- 
formed nuclei, following mainly the magic neutron num- 
bers N = 50 and 82. The sequence for the Gogny 
HFB calculation exhibits deviations from the other nu- 
clear models, visibly by the appearance of the nucleus 
fgCd with a large proton number and a jump to ^fRu 
in the sequence thereafter. The reason for this behavior 
is the comparatively low binding energy of the nuclei in 
this region, which is probably due to the calculation be- 
ing performed in the harmonic oscillator basis or to the 
parametrization D1S itself. 

The composition of the outer crust of neutron stars 
and the role of higher order corrections have been studied 
also in detail. The biggest influence of the higher order 
correction on the equation of state was found at lower 
densities, where the higher order corrections decrease the 
total pressure by a few percent. For larger densities the 
higher order corrections were almost negligible for the 
equation of state. Nevertheless, they led to small changes 
in the sequence of nuclei which appear at high densities 
> 10 11 g/cm 3 ). 

In both calculations no triaxially deformed nuclei were 
present in the sequence of nuclei. One can conclude that 
the effect of triaxial deformations is too weak to have 
significant consequences for the sequence of nuclei of the 
outer crust of neutron stars. 

In conclusion, the location of the drip line and the 
sequence of nuclei in outer crust of cold nonaccreting 
neutron stars seem to be rather robust predictions being 
nearly insensitive to the parameter set, the approxima- 
tion scheme and even to triaxial deformations for most 
of the state-of-the-art nuclear models. These findings are 
of importance for future rare isotope facilities as ISAC-H 
and FAIR at GSI Darmstadt, where binding energies of 
neutron-rich nuclei will be experimentally determined in 
the near future. 
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